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Abstract 

We report unusual jamming in driven ordered vortex flow in 2H-NbS2. Reinitiating movement 
in these jammed vortices with a higher driving force, and halting it thereafter once again with a 
reduction in drive, unfolds a critical behavior centered around the de-pinning threshold via 
divergences in the lifetimes of transient states, validating the predictions of a recent simulation 
study, which also pointed out a correspondence between plastic de-pinning in vortex matter and 
the notion of random organization proposed in the context of sheared colloids undergoing 
diffusive motion. 
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Djmamical transitions in driven non-equilibrium systems have been considered to exhibit 
signatures similar to equilibrium critical phase transitions, e.g., the divergence in length or time 
scales with characteristic critical exponent^'^'^'"^. Recent simulations and experiments on driven 
two-dimensional (2D)-conoidal systems have shown that the transition from a pinned state to a 
plastic flow state is associated with diverging time scales^, along with a characteristic bimodal 
velocity distribution^. The initiation of sliding motion in vortex arrays in type-II superconductors 
is analogous to the generic de-pinning transition in condensed matter, such as, de-pinning of 
charge density waves, colloids, Wigner crystals, magnetic bubbles, magnetic-domain walls, etc^. 
A widely investigated non-equilibrium driven situation in superconducting vortex matter is the 
plastic de-pinning of vortices^'*'^'^^'^^'^^'^^'^'^ wherein islands of pinned vortices coexist with 
filamentary channels of freely flowing vortices. 



Application of an external transport current normal to the magnetic field gives rise to a Lorentz 
force which de-pins the vortex state to reveal a variety of different driven phases, e.g., a moving 
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ordered Bragg glass state, the plastic flow state, etc. Usmg a 

new protocol, we have uncovered in a single crystal of a low Tc anisotropic superconductor 2H- 
NbS2, two Umiting ends of the de-pinning force, viz., de-pinning of an ordered pristine pinned 
vortex state (having a lower critical current density, J[ ) and that of a disordered vortex state 

(having a higher critical current density, j''^ ). In the present context, the disordered vortex state 
gets nucleated as the free flow (FF) in the driven ordered vortex state is jammed either by 
attempting to steadily accelerate or by waiting for long time. We also report on the discovery of 
vortex-velocity fluctuating between two extremes near the threshold limit of de-pinning the 
jammed, disordered vortex state. The life-time of the moving state in the fluctuating mode is 
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found to diverge on approaching from above as well as below it, validating the prediction of 
a simulation study^, which drew connection between plastic depinning and the phenomenon of 
random organization"*. Few recent experimental studies^^'^^ in periodically driven vortex matter 

TO 

pursuing reversible - irreversible transition in periodically driven systems have also explored 
the above mentioned connection articulated in Ref. 5. 

The isothermal electrical transport measurements have been performed at T = 2.5 K in the 
standard four probe geometry, with the dc magnetic field (H) applied parallel to the c-axis of the 
hexagonal single crystal of 2H-NbS2, and with the dc current (7) flowing in its basal (ab) plane. 
The single crystals of 2H-NbS2 were grown by the standard vapour transport technique. The 
crystal used has dimensions of 0.9 x 0.9 x 0.045 mm^, with Tc of 5.8 K and the residual 
resistivity ratio, R(300 K)/R(10 K), of 35. When current / sent thorough a superconductor 
exceeds the threshold de-pinning current (4), a voltage V appears. We use a criterion of 
appearance of V of 5 to determine Ic. The relationship V=Hud, where d is the separation 
between the voltage contacts (~ 0.3 mm in the present case) helps in estimating the mean 
velocity (m), with which the vortices drift when / > Ic, 

Figure 1 presents V-H curves with / = 15 mA at 2.5 K, for two different field ramp rates, 
Note first that the de-pinning with a given / (i.e., development of voltage) 

commences at a lower threshold field of 4 kOe for H -300 Oe/min, as compared to that near 13 
kOe for H = 100 Oe/min (corresponding to m = 23 cm/s) . The V-H response for H = 100 
Oe/min elucidates that Ic(H) is greater than the applied / of 15 mA for // < 13 kOe. After de- 
pinning at H = 13 kOe, the V-H gradually enhances to reach its normal state value at upper 
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critical field, Hc2. Considering that while ramping @ H =300 Oe/min, the de-pinning with 7=15 
mA commences at 4 kOe, the /c is > 15 mA only for H < A kOe. Thereafter, as the Ic{H) 

12 

progressively decreases further (below 15 mA), the voltage (proportional to a positive power 
of (7 - 7c(77)) continuously enhances from 4 kOe up to 10 kOe, after which the vortex flow is 
abruptly halted. This corresponds to a djmamic re-pinning, equivalent to a jamming transition in 
the driven vortex state. The V-H data in Fig. 1 establishes that the vortex states prepared with 
two different 77 have entirely different 7c(77) values for 4 kOe < 77 < 10 kOe. Above 10 kOe, the 
two V-77 curves in Fig. 1 nearly overlap, and de-pinning commencing at about 13 kOe and the 
voltage response thereafter is independent of 77 . 

31 * 

Figure 2 shows a set of V-7 curves for vortex states at 8 kOe prepared at 2.5 K via different 77 , 
ranging from 100 Oe/min to 250 Oe/min. These curves clearly identify two distinct 7c (at 77 = 8 
kOe) values, 7' and 7^ , as marked in Fig. 2. The inset (a) of Fig. 1 summarizes the variation of 
7c (77 = 8 kOe) with 77 . The higher Umiting value 7* (~ 56 mA) is obtained with 77 < 180 
Oe/min, while the lower limit 7' (~ 10 mA) identifies with 77 > 200 Oe/min. 

The V-7 curve for the state prepared with 77 = 250 Oe/min in Fig. 2 has the usual inflection 
feature/knee shape^^'^^ at 7 = 7'^, this crossover value^^ can be conveniently ascertained from 

dV/dl vs. 7 plot (cf. Fig. 14, Ref. [13]). Above , a quasi-linear V-I response sets in (implying 
steadily accelerating uniform flux flow (FF) regime), as 7 is gradually raised up to 88 mA. The 
state prepared with 77 = 200 Oe/min, also, has a low Id- 10 mA). However, in its quasi-linear 
flux flow regime above about 40 mA, it displays very large excursions, with V intermittently 
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dropping to nominal zero, corresponding to the jamming of vortex flow. The fluctuations persist, 
and nearly cease only for / > 75 mA. The presence of fluctuations at / > 60 mA can also be 
witnessed in the V-I curves for the states prepared with H = 150 Oe/min, where the initially 
prepared vortex matter has a higher critical current The higher crossover current limit, at 
which quasi-linear response sets in and large fluctuations in V largely cease, has been denoted 
by Icr (~ 83 mA) in Fig. 2. For H - 100 Oe/min, the V-I data are depicted for / ramped up to 
about 91 mA, and then back down to zero (indicated by the two oppositely directed short and 
thick arrows in Fig. 2). While ramping down, the V-I curve does not retrace its path below the 
crossover current /^^ • Larger excursions in V-I during ramp down commence as / is lowered 
below 56 mA {=Ic), however, they eventually cease, and the moving vortex state once again 
reorganizes into the disordered jammed vortex configuration, when 7-30 mA. On ramping up 
the current from zero value for the second time, the de-pinning commences once again only at 
the higher critical current /* , and not at the lower critical current, /[, . This observation contrasts 
the finding in weakly pinned 2H-NbSe2^'^ that a vortex state with the lower de-pinning critical 
current is nucleated upon halting a steady flowing state created out of driving a super-cooled 
(i.e., field cooled) disordered vortex state with a higher IdH). 

The inset (b) in Fig. 1 summarizes the variation of /^and with H at 2.5 K. Due to collective 
interaction effects between the vortices, the Ic(H) is known to depend inversely^"* on H. The field 
dependent lower I'dH) values are thus identified with the ordered vortex matter, and the nearly 
field independent (upto 10 kOe in the inset (b)) higher /* (H) values characterize the disordered 
counterparts ^'^"^'^^ corresponding to independent pinning of small bundles of vortices. 
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The large fluctuations sampled in the vortex-velocity in Fig. 2 motivated us to record the time 
series of V(t) response. Figure 3 collates V(t) plots atH=8 kOe (reached via H = 250 Oe/min, 
after initial zero field cooling at 2.5 K) for six current values. In Figs. 3(a) to 3(c), the initial 
situation on the average (at f = 0) is a moving vortex matter (with finite voltage level, above 5 

|aV) and the eventual condition after a long waiting time (j-'^) (ranging from -10^ to -10^ 
seconds) is a jammed state (i.e., V < 2.5 |aV). When / = 11.5 mA (cf. Fig. 3(a)), the notionally 
ordered vortex matter having a lower critical current (/[,- 10 mA) is just depinned, and the 
mean voltage level along with the average vortex velocity is small (~ 2 cm/s), but, after wait 
time^-;^ of about 1300 seconds, the V abruptly drops to nominal zero level (i.e., < 5 |aV). In Fig. 
3(b), at / = 28 mA (> crossover current ), the moving vortex state remains in quasi-linear 
flow mode, with small excursions anchored around <V> ~ 150 |aV, for about 3000 seconds. 
However, at / = 40 mA (cf. Fig. 3(c)), is only about 350 seconds. 

The ^'^ values for l[ <I < I'l are plotted in Fig. 4, and they display a non-monotonic behavior. 

The progressively faster moving vortex matter flow reorganizes into a disordered vortex 

configuration in a shorter time at / > 30 mA, with j'^ reducing to just about 50 seconds at / = 50 

mA. We believe that enhancement in the driving current fuels a competition between the 
generation of defects''^'^^'^^ i.e., dislocations, and their subsequent annealing out in the driven 
vortex state . One may therefore associate a net rate at which the topological defects, like 

dislocations, proliferate in the vortex state system. Therefore, ^^^is a measure of the drive (7) 
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dependent time scale needed for a threshold density of defects to build up in the moving lattice 
at a given field in a given sample^^ beyond which the stability of the flowing vortex state is 
significantly compromised (which in the present case occurs beyond 30 mA), and the flow 
spontaneously reorganizes into disordered jammed vortex configuration. The vortex- motion 
induced crossover to a higher Ic{H) state^^ reported earlier in YNi2B2C is notably different from 
the present result of a sudden change to a disordered state in 2H-NbS2, due to similar field 
dependences of the higher Ic{H) and lower IdH) states in YNi2B2C. These two field dependent 
Ic{H) states in YNi2B2C^^ are thus identified with ordered vortex lattices of different symmetry 
possible to observe in this system. 

The disordered configuration realized after a wait time of j'^ can be de-pinned by enhancing the 

current beyond /* (= 56 mA). Figure 3(d) shows the V{t) at / = 60 mA. Note first that initially 
(i.e., up to about 50 seconds), the voltage level rapidly fluctuates between 20 |aV and 60 |aV and, 
thereafter, the upper limit of fluctuations enhances to reach up to about 190 |aV. The 'state of 

fluctuations' lives up to about 1000 seconds (marked as ^^), and then suddenly a steady free 

flow (FF) (corresponding to V ~ 160 |aV in Fig. 3(d)) emerges. The measured for different / 

(>/c) are also plotted (as open circles) in Fig. 4. Note that values decrease rapidly as / 

progressively increases above . The identification of fluctuations and the divergence of j[ as 

/-^/* in Fig. 4 are two of our key observations. The notion of divergence stands vividly 

illustrated in the V{t) plot for / = 56 mA in Fig. 3(e). At / , the large fluctuations in V{t) 

range between just above zero level (closer to a pinned disordered state) to the high level of 
about 490 |aV. Figure 3(f) shows that vortex matter initially driven into a steady flow with 7=91 
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inA (>/c), when slowed down by decreasing the current to 43 mA (/ </c), maintains its 
flowing mode at the reduced current only for a short duration (-30 seconds). Thereafter, it 
fluctuates and eventually attains the disordered jammed vortex configuration, j-^il) measured 

for different / (<Ic), reduced from initial current of 91 mA, have also been plotted (as filled 
circles) in Fig. 4. 

Figure 4 shows that diverges from both below as well as above , viz., the back and forth 
transformation between a steady flowing and a jammed vortex state is through the intervening 
fluctuating states, whose transient life times diverge upon approaching /^from either side of 

it^"^. The inset in Fig. 4 shows a fit of the ^[(I) data to the relationship^ T ocl (/ - ) f , with ~ 
- 1.60 (+0.12) for the disordered vortex matter at 8 kOe at 2.5 K . The diverging timescales of 
large amplitude fluctuations (fl^if)) in close vicinity of suggest that, like few other non- 
equilibrium systems, the FF vortex state can undergo a random organization into a pinned 
disordered configuration, whose de-pirming attribute is like a djmamical phase transition in a 

2 3 4 5 25 

driven non-equilibrium system ' ' ' ' . Recent studies suggest the parameter fi, which is akin to 
a critical exponent^, could imbibe the notion of universality class for dynamic transitions in 
diverse types of driven systems. In a different context^^'^^, a case was made for the plastic 
depinning as a djmamical transition in 2H-NbSe2 (see Fig. 25 in [13]). However, the present 
results of characteristic diverging time scales on both sides of the de-jamming threshold in 2H- 
NbS2 establish the correspondence between the plastic depinning in the driven vortex state and a 
(djmamical) random organization occurring between a non-fluctuating (pirmed or immobile) and 
a fluctuating (flowing) state in driven (diffusive) colloidal systems"^'^. 
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The bi-modal character of the moving phase fluctuating between two hmiting values of the 
vortex velocities also echoes the behavior reported in the plastic flow regime of driven colloidal 
crystals'", and in simulation studies'" on driven vortex matter. Analysis of the probability 
distributions of vortex-velocities observed by us^^ shows that a fraction of the slower moving 
vortices coexist with a faster moving fraction. The analogue of vortex phase fluctuating between 
two extreme values of vortex velocities was not identified in the context of colloidal matter^, 
which displayed the divergence characteristic in the transient times. The observation of negative 
values of vortex-velocities, opposite to the direction of drive (see, for instance, Ref. 35 and cf. 
Fig. 3) in the disordered jammed vortex state is also a new observation and justifies the jammed' 
nomenclature assigned to the disordered state in 2H-NbS2. Similar negative (entropy 
consuming) events have also been observed and analyzed through the non-equilibrium steady 
state fluctuation relation, in the sheared jammed state in surfactant-based hexagonal phase of 
cylindrical micelles^^. 
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funding support from Department of Science and Technology of Govt, of India. 
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Fig. 1 . Gorky Shaw et al, 2011. 

Fig. 1: (Color online) Voltage drop versus field, V-H, across the single crystal of 2H-NbS2 for 
// II c and at 2.5 K for / = 15 mA with H = \00 Oe/min (open circles) and H = 300 Oe/min 
(filled circles), respectively. The two arrows identify the threshold field values (13 kOe and 4 
kOe, respectively) at which the V{H) rises above 5 |^V. The V{H) for H = 300 Oe/min displays a 
sudden transition from a high voltage (i.e., moving state) to the nominal zero value (i.e., jammed 
state) at // = 10 kOe (« Hc2). The inset (a) is a plot of 4 vs. dH/dt at 8 kOe and the inset (b) 
displays variation of and vs. // (II c) at 2.5 K. 
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Fig. 2: (Color online) The V-/ responses in 2H-NbS2 for some selected vortex states at 8 kOe {H 
II c), prepared via different// , ranging from 100 Oe/min to 250 Oe/min. The critical currents, 
![ and /^' , stand identified. The limiting values and /*. , above which V-I curves are quasi - 
linear are also marked. For vortex state generated via // = 100 Oe/min, a reverse leg of the V-I 
curve has also been shown. The cross marks (x) on the V-I curve for H = 250 Oe/min identify 
the three current values at which the time series measurements on the voltage are presented in 
Figs. 3(a) to 3(c). 
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Fig. 3 : (Color online) The panels (a) to (f) display time series measurements (data points shown 
as red open circles with the symbols joined by blue line) of voltage drop (cc vortex velocity) at 
different / in a single crystal of 2H-NbS2 at H = S kOe (II c) and T = 2.5 K, (a) / = 11.5 mA 
(>/^' ), (b) / = 28 mA (>/^,) and (c) / = 40 mA, (d) / = 60 mA, i.e., just after de-pinning the 

jammed/disordered state, (e) / = 56 mA (~ /'' ), i.e., at the threshold of de-pinning the disordered 
state, and (f) / = 43 mA, i.e., on reduction of current from free flow (FF) state at 1= 91 mA to / = 
43 mA. The moving state (<V> is greater than 5 |_iV) can be seen to transform into the jammed 

state {<V> ~ 2.5 |iV) after a lifetime of j-'^in panels (a) to (c). In panels (d) and (f) the 

fluctuating state can be seen to settle down to the free flow mode after a transient time of 2", • In 

panel (e), the fluctuating state seems to persist forever {we recorded data up to about 4100 
seconds). Note also the overshoot of the vortex velocities below zero, i.e., to negative values in 
different panels. 
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Fig. 4: (Color online) Variation of 2",' and 2",^ with current / for vortex states at 8 kOe (lie) and 

T= 2.5 K. Three encircled data points on the j-'^ (I) plot represent the typical / values at which 
the time series data have been displayed in Figs. 3(a) to 3(c). An inset panel shows plots of log 
-3-,' vs. log 1(7- /c)l for two sets of data, for (i) I > /'' (open circles) and (ii) I < it (filled 
circles), respectively. The straight line drawn amounts to a power law relationship, 
tI (^-^c) , with exponent, -1.60 ± 0.12. 
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